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15 Abstract 

1 6 Both sensor calibration and satellite retrieval algorithm play an important role in the ability to 

17 determine accurately long-term trends from satellite data. Owing to the unprecedented 

18 accuracy and long-term stability of its radiometric calibration, the SeaWiFS measurements 

19 exhibit minimal uncertainty with respect to sensor calibration. In this study, we take 

20 advantage of this well-calibrated set of measurements by applying a newly-developed aerosol 

21 optical depth (AOD) retrieval algorithm over land and ocean to investigate the distribution of 

22 AOD, and to identify emerging patterns and trends in global and regional aerosol loading 

23 during its 13-year mission. Our results indicate that the averaged AOD trend over global 

24 ocean is weakly positive from 1998 to 2010 and comparable to that observed by MODIS but 

25 opposite in sign to that observed by AVHRR during overlapping years. On a smaller scale, 

26 different trends are found for different regions. For example, large upward trends are found 

27 over the Arabian Peninsula that indicate a strengthening of the seasonal cycle of dust emission 

28 and transport processes over the whole region as well as over downwind oceanic regions. In 
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contrast, a negative-neutral tendency is observed over the desert/arid Saharan region as well 
as in the associated dust outflow over the north Atlantic. Additionally, we found decreasing 
trends over the eastern US and Europe, and increasing trends over countries such as China 
and India that are experiencing rapid economic development. In general, these results are 
consistent with those derived from ground-based AERONET measurements. 


1 Introduction 

The impact of natural and anthropogenic sources of air pollution on Earth’s weather and 
climate systems and their long-term tendencies have gained increasing attention from the 
scientific community in recent years. Indeed, tropospheric aerosols not only perturb the 
radiative energy balance by interacting with solar and terrestrial radiation (Ramanathan et al., 
2001) but also by changing cloud properties and lifetimes (Rosenfeld et al., 2008). 
Furthermore, the Intergovernmental Panel on Climate Change (IPCC, 2007) report indicates 
that the aerosol cooling effect could have partially counteracted warming from greenhouse 
gas increases over the past few decades. Elowever, there are large uncertainties in the 
estimation of climate forcing from aerosols due to their complex nature and short lifetime. In 
order to achieve a better understanding of the spatial and temporal variability of aerosol 
distributions on both a regional and global scale, long-term satellite measurements of high 
fidelity are required. From this increased understanding we can then accurately evaluate 
aerosol effects in the climate models. 

Despite its importance, deriving the small signal that represents the long-term trend in aerosol 
distribution from the large seasonal cycles and natural year-to-year variability which result 
from large-scale meteorology is highly challenging. In order to use satellite measurements to 
determine aerosol trends with the required level of robustness, several key elements must be 
addressed. They include: 1) accuracy, precision, and stability of the sensor calibration, 2) 
potential for sampling bias due to insufficient temporal or spatial coverage, 3) potential for 
algorithm bias (i.e., differences in the cloud screening scheme, temporal changes in the 
aerosol or surface properties that deviate from assumptions used in the retrieval algorithm), 4) 
sufficient length of the data record, and 5) potential effects of any measurement time drift due 
to the diurnal cycle of atmospheric aerosol loading. Early satellite sensors, such as the 
Advanced Very FTigh Resolution Radiometer (AVHRR), provide a longer data record for 
trend detemrination compared to later, more advanced satellite sensors such as the Sea- 
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viewing Wide Field-of-view Sensor (SeaWiFS), Moderate Resolution Imaging 
Spectroradiometer (MODIS), Multiangle Imaging Spectroradiometer (MISR), and others. 
Flowever, due to the limited number of spectral bands and the lack of onboard calibration for 
these early sensors, the uncertainty associated with trends calculated using measurements 
from them can be significant. 

Among the suite of NASA’s Earth Observing System (EOS)-era satellites, the SeaWiFS 
instrument was primarily designed to measure ocean color; since its launch in 1997, it has 
been the major source for providing a comprehensive global data set of such measurements. 
Flowever, the unprecedented set of well-calibrated radiances measured in the wavelength 
range from the visible (412 nm) to the near infrared (865 nm) also make it well-suited to 
providing information about atmospheric aerosols [McClain et al., 1998], In this study, we 
will attempt to estimate trends in aerosol optical depth (AOD) at 550 nm over land and ocean 
using SeaWiFS measurements from 1997 to 2010 using a new SeaWiFS aerosol dataset that 
has been produced as part of NASA’s MEaSUREs project. 

We begin in Section 2 by evaluating the performance of the sensor calibration and aerosol 
retrieval algorithm throughout the entire SeaWiFS mission. Section 3 describes our 
methodology for trend determination. This is followed, in Section 4, by a discussion of the 
resulting annual and seasonal trends from SeaWiFS and comparisons to ground-based 
measurements from the Aerosol Robotic Network (AERONET, Holben et al., 1998). A 
summary and conclusions based on our study are given in the closing section. 


2 Satellite data and their uncertainties for trend analysis 

Since 90% of the top-of-the-atmosphere (TO A) radiance received by a satellite sensor over 
ocean comes from contributions of the atmosphere such as Rayleigh and aerosol scattering, 
less than 10% of the signal is attributable to the ocean itself. Consequently, the level of 
accuracy, precision, and long-term stability needed in radiometric measurements from such a 
sensor is much higher for ocean color retrievals than those needed for corresponding land or 
atmospheric products. As a result of rigorous on-orbit calibration techniques developed using 
a combination of lunar and solar diffuser measurements and vicarious procedures, the NASA 
Ocean Biology Processing Group (OBPG) Calibration and Validation Team achieved an 
extraordinary level of absolute accuracy (1-2%), precision (0.1%) and, most importantly, 
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long-term stability (0.3% over the 13-year mission) for the TOA radiances from SeaWiFS 
(Eplee et al., 201 1). 

The SeaStar platform carrying SeaWiFS was launched with a nominal center-of-swath 
equatorial overpass time of around noon for the daytime (descending) node. This drifted 
slowly through the mission into a later orbit, although remained before 12:30 pm until 2005; 
the time reached a maximum around 2:40 pm in July 2010, after which an orbit raising 
maneuver was performed to return the orbit slowly closer to noon. Smirnov et al. (2002) 
investigated the diurnal cycle of AOD at various AERONET sites. Over this range of local 
times, the change in mid-visible AOD was found to be of order of a few percent or less of the 
daily mean AOD for most aerosol types; the most rapid diurnal variability was generally in 
early morning or late afternoon. This suggests that the drift in SeaWiFS orbit time is unlikely 
to introduce significant artifacts into derived trends. 

The SeaWiFS product used in this study is based upon the SeaWiFS v003 dataset that is 
available from http://disc.gsfc.nasa.gov/dust. This dataset spans the time period September 
1997 to December 2010 and is the result of an effort funded by NASA’s MEaSUREs project 
to produce a long-term aerosol data record from SeaWiFS. This aerosol retrieval system 
utilizes a combination of the Deep Blue algorithm over land (Hsu et al., 2004, 2006) and the 
SeaWiFS Ocean Aerosol Retrieval (SOAR) algorithm over ocean (Sayer et al., 2012a). The 
Level 2 products are derived with a horizontal resolution of approximately 13.5 km at nadir 
and include spectral AOD, Angstrom exponent, and single scattering albedo (for dust 
aerosols) over land and spectral AOD, Angstrom exponent, and fine mode fractional aerosol 
volume over ocean. Level 3 (gridded) products are available for daily and monthly averages at 
both 0.5° x 0.5° and 1° x 1° resolutions. Figure 1 shows an example of the resulting 
seasonally-averaged AOD at 550 nm throughout the 13 year mission using this newly- 
developed product. In this study, we used the Level 3 1° x 1° spatial resolution monthly mean 
SeaWiFS AOD data at 550 nm to calculate aerosol trends that could be efficiently compared 
with other satellite products. 

A detailed validation of the global SeaWiFS AOD using the collocated v003 Level 2 products 
against the cloud-screened and quality-assured Level 2 AERONET measurements (Holben et 
al., 1998, Smirnov et al., 2000) and other satellite products is presented for land in Sayer et al. 
(2012a) and for ocean in Sayer et al. (2012b). Briefly, on a global basis, over one standard 
deviation (68%) of AOD matchups were found to agree within an absolute expected error of 
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0.03+15% over ocean and 0.05+20% over land at 550 nm. Since a change in the quality of the 
retrieval with time could mask or amplify any true change in AOD, it was also imperative to 
make an assessment of the temporal stability for this dataset. We therefore calculated the 
differences in 550 nm AOD between SeaWiFS and AERONET data as a function of year over 
ocean and land separately; the results are depicted in Figure 2 (methodology and general 
validation results are presented by Sayer et al., 20 12a, b). Only the best quality SeaWiFS data, 
with a quality assurance (QA) of 3 over land and a QA of 2 and 3 over ocean, are included in 
our analyses (as recommended in the aforementioned references). Also, in order to ensure the 
consistency of AERONET measurement locations throughout the SeaWiFS lifetime, only 
AERONET sites with records lasting for most or all of the SeaWiFS time series are 
considered, as otherwise an apparent change in comparison statistics could simply be the 
result of changes in regional sampling. 

For retrievals over land, the stations used in our analysis are Alta Floresta, Avignon, 
Banizoumbou, Beijing, Bondville, Bratts Lake, the CART site, Dakar, Dalanzadgad, GSFC, 
HJ Andrews, Ilorin, IMS METU-Erdemli, Ispra, Kanpur, MD Science Center, Mongu, Sede 
Boker, Sevilleta, Skukuza, Solar Village, and Wallops. Figure 2 shows that the median and 
mean AOD difference between SeaWiFS and AERONET for these sites is generally small 
and negative; therefore, for a typical case, SeaWiFS is likely to have a slightly low bias. The 
lines are reasonably stable over time, suggesting that the quality of the dataset does not 
change throughout the SeaWiFS mission. The intercept of the best-fit line of the mean error 
against time is -0.016 (±0.007), with a gradient of 0.0003 (±0.0009) per year (both in units of 
absolute AOD). Therefore, the temporal bias is not statistically different from zero at the 
95% level (quoted fit uncertainties are the one-sigma, i.e. 68%, confidence limits). Variation 
from year-to-year appears smooth, with no abrupt changes. 

Over ocean, fewer sites are available, particularly for the early part of the time series, so only 
the years 1999-2010 are considered, from the AERONET sites at Arica, Capo Verde, Dakar, 
Flelgoland, IMS METU-Erdemli, La Parguera, Nauru, Shirahama, and Venise. The intercept 
of -0.008 (±0.005) and gradient of -0.0014 (±0.0008) per year of the best-fit line are not 
significantly different from zero at the 95% confidence level (the gradient is significant at the 
90% level, but is small in magnitude). 

These results give confidence that the SeaWiFS aerosol dataset is suitable for the analysis of 
changes in AOD over this period. This does, of course, rely on the assumption that the 
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evaluation at these sites is a good representation of the whole globe, although in the absence 
of long-term aerosol measurements at every location, this cannot be assessed directly. 

3 Methodology 

In order to show changes in spatial and temporal distribution of aerosol loading, linear trends 
were calculated based on the 1° x 1° spatial resolution monthly mean SeaWiFS AOD data, 
from September 1997 to December 2010. These monthly datasets are calculated as the simple 
mean of the daily mean AOD for each month; the daily mean is computed from all those 
retrievals of sufficiently high QA (2 or 3 over ocean, 3 over land). In order to remove the 
large influence of the annual cycle, the data were first deseasonalized by obtaining the time 
series of AOT anomaly for each grid cell. The anomaly is defined as difference between 
monthly variation and climatology (i.e. the month-wise average over all the years). A least- 
squares fit was applied to the AOD anomaly time series to get the slope coefficient of the 
linear regression that represents annual trend (AOD/year). 

The statistical significance of the linear regression was assessed using the analysis of variance 
procedure which involves testing the F statistic and associated p value. The F statistic is 
referred to as the ratio between mean-squared error of regression and mean-squared error of 
residuals. The associated p-values (less than 0.05) are indicated as "dots" on global trend 
distribution maps for regions that are statistically significant at the 95% confidence level, and 
the null hypothesis (that the regression coefficients are not statistically different from zero) 
can be rejected. 

As least-squares regression can be sensitive to outliers, autocorrelation, and heteroscedasticity 
of errors in the input data, linear AOD trends were also estimated using the non-parametric 
Theil-Sen (aka ‘median gradients’) method, which is an estimator robust to these difficulties 
(Theil, 1950, Sen, 1968). Global trends calculated using the Theil-Sen method were within 
approximately 10% of the least-squares regression values for all the sensors considered, 
suggesting the least-squares model is appropriate for this situation. As the least-squares 
regression technique more readily provides uncertainty estimates, only the least-squares 
results are used hereafter. 
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4 Results and discussions 


4.1 Correlations between SeaWiFS AOD anomaly with ENSO and NAO 

In order to estimate global and regional trends in AOD, we need to first examine the inter- 
annual variability of the aerosol distribution and its relationship with large-scale climatic 
variations such as the El Nino/La Nina-Southern Oscillation (ENSO) and the North Atlantic 
Oscillation (NAO). These large-scale meteorological processes could perturb the aerosol 
distributions through changes in aerosol emissions by altering the precipitation and soil 
moisture or in aerosol transport via varying the atmospheric circulation. The correlation 
coefficients (R) are calculated between deseasonalized monthly SeaWiFS AOD anomalies 
and the Multivariate ENSO Index (MEI; Wolter and Timlin, 1993, 1998) from 1997 to 2010. 

As shown in Figure 3(a), the most prominent regions of high correlation with ENSO 
signatures are in the tropics near Indonesia and Central America. In Indonesia/Bomeo, fires 
are primarily used to clear vegetation waste associated with deforestation for agricultural 
expansion and usually peak in the fall and spring season, when the air is relatively dry and 
precipitation is low. These smoke plumes, which can often stretch over hundreds of 
kilometers, are readily observed by satellites. The times series of SeaWiFS monthly AOD 
anomaly along with ENSO index is shown in Figure 4 and indicates that the intensity of 
resulting smoke is predominantly linked to El-Nino-induced droughts over this region; strong 
peaks are seen in the years 1997-1998, 2002, and 2006. 

The amount of precipitation affects the interannual variability of biomass burning activities 
and dust emission through changes in surface moisture; since these effects could be 
accumulative, we also calculate the correlation coefficients between SeaWiFS AOD anomaly 
and the ENSO index with various time lags. Figure 3b and 3c show these AOD anomaly 
correlations with ENSO index leading by six months and one year, respectively. We note that 
the intensity of biomass burning smoke in Central America is sensitive not only to the ENSO 
induced rainfall changes during the peak burning season, but also to those accumulated 
throughout several months prior to the start of burning (i.e., on a 6-month time lag as shown 
in Fig. 3b). The lag time for South America appears to be even longer (i.e., 6 months to a 
year) compared to that for Central America. Such time lags between ENSO and fire season 
severity in South America have also been reported by Chen et al. (2011), with significant 
lagged correlations emerging on a 3-7 month time scale. However, there seems to be little 
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time lag observed over Indonesia between AOD anomaly and ENSO-induced drought. This is 
probably indicative of the differences in agriculture/land clearing practices between South 
America, Central America, and Indonesia. 

The results shown in Figure 5 indicate that summertime Saharan dust loadings highly 
correlate (R > 0.5 at 95% confidence) with the ENSO index from 6 months earlier over dust 
source regions in northern Sahel, as well as over the Atlantic along the dust transport pathway 
reaching as far as Barbados, West Indies. This is consistent with the findings of Prospero and 
Lamb (2003) and Chiapello et al. (2005), who reported the close link between dust 
concentrations observed in Barbados followed by the occurrence of major ENSO events and 
the associated drought in Sahel. In contrast to the Sahara, which receives little rainfall each 
year, dust emissions over the Sahel are apparently more influenced by the year-to-year 
variability of the precipitation, as expected. Besides the transport route from the Sahel to 
Barbados, the high ENSO/AOD correlation results shown in Figure 5 also indicate the 
possibility of a widespread association between ENSO and dust transport in the western and 
northern parts of the North Atlantic along the edge of the Saharan dust outflow regions. This 
anomalous pattern also coincides with the main development region (MDR), primarily 
between 10°-20°N, of tropical cyclone activity in the Atlantic. This is likely to be linked to 
the effects of ENSO induced changes of the atmospheric dynamics and circulation on the 
summer time Saharan dust transport into the Atlantic. 

The effects of the NAO on the Saharan dust distribution are also investigated as shown in 
Figure 6. The NAO is one of the dominant factors of winter climate variability in the North 
Atlantic region, spanning central North America and Europe. The NAO winter index used in 
this study is based upon the difference of normalized sea level pressure between Lisbon, 
Portugal and Stykkisholmur, Iceland (Hurrell, 1995). Based upon the results in Figure 6, a 
sizeable and statistically significant (at 95% confidence) correlation (R>0.5) between is 
apparent between the SeaWiFS AOD anomaly and NAO index over extensive areas 
encompassing the western part of the North African continent into the North Atlantic. In 
addition, we also observe a secondary branch of positive correlation that is statistically 
significant along the pollution outflow of the eastern US across the most northern edge of the 
Atlantic. We found no statistically significant correlation over this region for seasons other 
than winter. This is consistent with the findings of Ginoux et al. (2004) and Chiapello et al. 
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(2005) based upon the 17-year record of TOMS data (1979-2000, except for 1982 and 1 993— 
1996). 

These results suggest that, while the NAO has a large influence on the wintertime dust 
distribution over the Atlantic, significant variability of the AOD anomaly in this region is 
controlled by the ENSO during the summer months. We also note that large scale 
meteorological events such as ENSO and NAO appear to have larger impacts on the 
interannual variability of aerosol distribution across the tropics, either in the dust-laden or 
biomass burning smoke-dominated regions, than in the extratropical zones. Therefore, trend 
determination in the tropical regions is more susceptible to the influence of these large-scale 
climatic driving forces and consequently more difficult to derive, as compared to those 
outside the tropics. 

4.2 Global trend 

Figure 7 shows comparisons of the time series of absolute AOD and AOD anomaly averaged 
over global ocean from SeaWiFS with those from MODIS Terra (data record beginning in 
2000) and Aqua (beginning 2002). Although there is an offset between MODIS Terra and 
Aqua in terms of absolute values of AOD, the short-term tendency during the overlapping 
period is similar between the two (Zhang and Reid, 2010). This comparison shows that the 
interannual variability of the SeaWiFS AOD anomaly follows both MODIS Terra and Aqua 
until mid- 2009. From 2009 onwards, the SeaWiFS AOD anomaly is higher than that from 
MODIS, resulting in a more positive trend over global ocean compared to the MODIS values. 
According to Figure 2, the comparisons of SeaWiFS retrieved AOD at 550 nm with the 
AERONET values do not indicate any drift from the zero line after the mid 2009. 

It is worth pointing out that long-term global ocean AOD measurements based on AVHRR 
data have been shown to have a downward tendency for the period 1982 - 2005 [Mischenko 
et al., 2007; Zhao et al., 2008]. Although revised calibration procedures applied to AVF1RR 
radiances have broadly resulted in comparable AOD when compared to MODIS [Zhao et al., 
2008; Fi et al., 2009], the short-term trends of AOD in the period of overlap between MODIS 
and AVHRR are opposite in sign (positive for MODIS and negative for AVHRR) for the 
period 2000 - 2005 (Fi et al., 2009; Zhang and Reid, 2010). This downward trend in AVHRR 
data continues until 2009, based on our analysis of the PATMOS-x AOD data (not shown), 
which is in contrast to the trends observed from SeaWiFS and MODIS from this study. 
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The difference in sign associated with trends in the past decade could partly be due to the 
overall calibration accuracy/orbital drift in the AVHRR data record. Thomas et al. (2010) 
compared changes in AOD between Along-Track Scamiing Radiometer 2 (ATSR-2) and 
AVHRR data for 1995-2001. They found that the global downward trend in the AVHRR 
AOD was largely driven by the southern oceans, for which the latitudinal range of sampling 
during daylight hours changed significantly during the study period, due to drifts in the 
AVHRR platform’s orbit. Even for those regions where daylight availability is not an issue, 
diurnal variability in cloud cover can lead to changes in the availability of clear-sky pixels for 
AOD retrieval as the satellites drift, resulting in further sampling bias (e.g. Devasthale et al., 
2012). The drifts in orbit time for the AVHRR sensors were generally larger than for 
SeaWiFS, and when compositing multiple AVHRR records together to create a longer time 
series, there are time discontinuities when one satellite is replaced by the next. Additionally, 
other factors contributing to this noted difference may be related to data sampling issues 
arising from different cloud-screening schemes employed by AVHRR and SeaWiFS/MODIS. 

Also shown in Figure 7 is the time series of over-ocean AOD and anomaly from the data 
assimilation (DA) grade MODIS-Terra dataset, which is an attempt to create a product 
suitable for DA through bias correction and additional quality screening of the MODIS 
Collection 5 aerosol product (Zhang & Reid, 2006; Shi et al., 2011). Only the period 2000- 
2009 is available. These time series track the SeaWiFS record very closely. This is consistent 
with Sayer et al (2012), who found, for collocated daily DA MODIS-Terra and SeaWiFS data 
over ocean, a correlation of 0.87 and bias of -0.006 between the two (SeaWiFS lower). In 
contrast, the Collection 5 MODIS Terra had a slightly lower correlation with SeaWiFS (0.80), 
and SeaWiFS a more negative bias (-0.027). Only Terra data are shown; the offset between 
the two MODIS sensors becomes smaller in the DA-MODIS product (Shi et al., 2011; Sayer 
et al., 2012). 

Zhang & Reid (2010) found, using a similar methodology to this study, an annual trend of 
0.001 yr" 1 in global-average over-ocean DA MODIS Terra AOD (and a similar figure noted 
for Aqua). This is in good agreement with the results of this study from both the DA MODIS 
and SeaWiFS records. However, after performing a correction to both the Terra and Aqua 
data to account for a drift in the quality of the validation statistics against AERONET since 
2005 (Fevy et al., 2010; Zhang & Reid, 2010), the magnitude of this trend decreased by a 
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factor of three, and it became statistically insignificant. Their best estimate is therefore that 
there has been no significant change in global-average over-ocean AOD during this period. 

The resulting annual trends of SeaWiFS AOD anomaly over both land and ocean are derived 
for each 1° x 1° grid cell and are shown in Figure 8. Points marked with dots indicate trends 
significant at the 95% confidence level (p-value < 0.05). For comparison, we also calculate 
trends based upon AERONET measurements over selected sites using the same methodology, 
where sufficient and more-or-less continuous data records were acquired over a similar time 
period to SeaWiFS. These AERONET values and the corresponding SeaWiFS trends are 
listed in Table 1. Also, to help visualize the SeaWiFS and AERONET trend comparisons in 
the context of the surrounding regions and over the whole globe, we superimpose the 
locations of AERONET sites on top of the SeaWiFS annual trend map in Figure 7. The black 
and pink boxes and symbol (x) represent negative, positive, and negligible trends from 
AERONET data, respectively. 

In general, the annual trends from 1998 to 2010 are small over most of the world. However, 
for the mineral dust dominant regions, strong positive trends are seen over the Arabian 
Peninsula, while negligible trends of no statistical significance exist over most of the desert 
regions in China. In contrast, there is a negative tendency in the surrounding regions of the 
western part of North Africa as well as the Atlantic. This downward trend, which is most 
likely associated with the decreasing export of Saharan dust from North Africa into the 
Atlantic, is also reflected in the AERONET data at Capo Verde. For fine mode anthropogenic 
aerosols, positive trends are observed in the fast developing countries over India and 
eastem/southem China, most likely associated with an increase in urban/industrial pollution. 
Negative trends are also seen over eastern US and Europe, which are consistent with the 
AERONET trends at GSFC/COVE, and IFT-Leipzig, Germany, respectively. The downward 
trends over Europe and the eastern US may partly be due to a decrease in pollution aerosols 
associated with emission control. For biomass burning regions, although there seems to be 
considerable downward trend over the Central America, the estimated trends are not 
statistically significant for other areas such as South America, North America, Africa, and 
Southeast Asia. 

Finally, it is evident from the results shown in Table 1 that, except for Capo Verde and 
COVE, most of the island/coastal AERONET measurements exhibit either positive or 
negligible trends over the tropical Pacific (Tahiti and Nauru), South Atlantic (Ascension 
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Island), the Caribbean (La Parguera), and the northwest Pacific (Shirahama). The 
corresponding trend values derived from the SeaWiFS records seem to be comparable to 
those from the AERONET data. 

4.3 Seasonal and regional trends 

Since, as shown in Figure 1, there is strong seasonality in both the emission intensities of 
natural and anthropogenic aerosols over various source regions and their impact on large- 
scale aerosol distributions, trend analyses were performed separately for each season. These 
seasonal trends were calculated over each 10x10 grid for the DJF, MAM, JJA and SON time 
spans and are presented in Figure 9. The corresponding time series of seasonal mean AOD 
anomaly over northern India (20° - 30° N, 75° - 85° E), eastern China (30° - 40° N, 110° - 
120° E), the Arabian Peninsula (10° - 35° N, 35° - 60° E), Europe (43° - 55° N, 0° - 30° E), 
and the eastern United States (30° - 45° N, 70° - 90° W) are shown in Figure 10. In order to 
have a reasonable temporal distribution in the analysis, trends were computed based on the 
time series of individual months for a particular season (stacked together), instead of 
calculating a seasonal mean trend. In other words, a slope coefficient was estimated from a 
time series consisting of DJF 1997, followed by DJF 1998 and so forth, thus from a total of 39 
data points. Linear trends and their significance are shown in Figure 9 in the same fashion as 
in Figure 8. 

It is apparent from the results shown in Figure 9 that there are strong positive trends over both 
the Arabian Peninsula as well as the surrounding Arabian Sea, especially during spring and 
summer seasons. Strong upward trends are also found over the adjoining Persian Gulf. Thus, 
it is likely that the increasing aerosol load over the Arabian Peninsula, together with the 
increase over the northern Arabian Sea, has contributed to upward trends in AOD over 
surrounding oceanic regions largely associated with the spring- and summer-time dust 
outflow. We also found supporting evidence of a systematically increasing aerosol loading 
tendency from AERONET observations at Solar Village in Saudi Arabia, as well as a 
decreasing trend in Angstrom exponent (440 nm - 870nm) for the period 1999-2010 (Figure 
11). Recently, Yoon et al. (2011a) also reported an increasing (decreasing) AOD (Angstrom 
exponent) for the period 2001-2008, based on analysis of AERONET data from the Solar 
Village site. The upward trend in AOD from AERONET data is 0.018 yr" 1 for the 11 -year 
period shown in Figure 11, while a positive trend is also obtained from the SeaWiFS data for 
the overlapping period over a l°xl° region collocated with the Solar Village site (Table 1). 
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Furthermore, we discern from the top panel of Figure 1 1 that the seasonal cycle of aerosol 
loading has strengthened on an approximate inter-annual scale, with measurable decrease in 
the Angstrom Exponent. These observations lead us to posit that the seasonal cycle of dust 
emissions (loading as well as outflow) has strengthened in the past 11-13 years spanning the 
AERONET and SeaWiFS data record over this region. The region-wide trend from the 
SeaWiFS AOD data record (September 1997 to December 2010) is estimated to be 0.009 yr" 1 
(±0.00076) when averaged over the Arabian Peninsula (10°N - 35°N, 35°E - 60°E). As clearly 
noted in the spatial trend distribution (Figure 8), this large regional trend is an order of 
magnitude higher than the global mean trend (land and ocean) as well as the global mean 
trend only over land. 

The global mean trends are weakly positive with large variability and are only representative 
for the relatively short time scale of the data record (a little over 1 3 years) and therefore are 
deemed to be negligible or not significant. However, it is important to emphasize the large 
sustained upward trends observed over the Arabian Peninsula and the need to monitor future 
trends, and/or construct a longer-term aerosol record in retrospect to better quantify the 
changes in aerosol loading, causes and regional (and downwind) climate impacts. 

In contrast to the strong upward trends over the Arabian Peninsula, there appears to be 
negative trends over the Saharan dust dominant regions in western African and the North 
Atlantic, particularly during the winter season. However, these patterns in the trends seem to 
coincide with the AOD-NAO correlation map depicted in Figure 6. Therefore, the negative 
trend in the SeaWiFS AOD anomaly during the winter months is likely a reflection of the 
downward trend in the NAO strength from 1997 to 2010. Similarly, the locations of negative 
trends in the summer AOD anomaly over the Sahel and across the Atlantic also mimic the 
AOD-ENSO correlation shown in Figure 5. Consequently, the large-scale meteorological 
processes such as NAO and ENSO largely regulate the export of the Saharan dust in the 
Atlantic. Overall, trends observed over the North Africa exhibit relatively smaller trends and 
weak seasonal dependency. However, there seems to be a significant positive trend over the 
Sahel region during springtime which is also reflected in the AERONET data at 
Banizoumbou. 

Over East Asia, the aerosol loading appears to exhibit negligible trends of no statistical 
significance in eastern China during the active dust outbreak season (i.e., springtime); the 
corresponding AOD anomaly time series is depicted in Figure 10 (top-right panel). The 
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amount of mineral dust reaching this region seems to decrease from 1998 to 2004-2005, but 
became quite high again in the later years (2006-2010). Overall, the insignificant trends for 
the dust-dominant season over eastern China are due to large inter-amiual variability in the 
observed aerosol loading. In South Asia, the regions over northern India along the Indo- 
Gangetic plain are often influenced during the pre-monsoon season with heavy dust plumes 
transported from the Thar Desert in northwestern India (Gautam et al., 2011). Based upon 
Figure 10 (top-left panel), these dust activities seem to increase significantly from 1998 and 
peak around 2003. However, since 2003 the dust intensity is seen to steadily decrease during 
the spring (MAM). Recent studies based on satellite observations from MISR and MODIS 
have also alluded to this decrease in dust loading over northern India from 2000 onwards 
(Dey and Di Girolamo, 2011; Kaskaoutis et al., 2011). However, since the SeaWiFS data 
record spans a longer time, our results indicate an overall slightly positive trend from 1998 to 
2010 over northern India. For the rest of the world, there are no statistically-significant trends 
over arid and semi-arid regions, including over Australia and the southwest US. 

For biomass burning regions such as Indonesia, South America, and in particular Central 
America, the aerosol trends are observed to be negative from 1998 to 2010. However, the 
year-to-year variabilities of smoke emission over these regions are most likely driven by the 
ENSO induced effects discussed in the previous section. Our calculations also show that there 
are no observed aerosol trends of statistical significance over Southeast Asia (i.e. Indo-China) 
as well as regions in North America, Europe, and Siberia where forest fires frequently occur. 

Over regions where anthropogenic fine mode aerosols are prevalent, our study shows that 
there are downward trends in AOD over the eastern US and Europe during the summer 
months (i.e. the peak season of air pollution build-up in these regions). However, for the 
rapidly developing countries in Asia, the AOD tendencies are positive. In particular, during 
the heavy winter haze months, significant upward trends are observed in the regions of 
eastern and southern China, the East China Sea, and in northern and central India and the 
surrounding oceanic regions, particularly over the Bay of Bengal. In contrast to the weak 
tendencies in the boreal spring/pre-monsoon months that are largely due to inter-annual 
variability of dust transport, the relatively stable conditions during dry winter/post-monsoon 
period favors accumulation of aerosols/pollutants. Hence, the upward trends are most likely 
due to increase in regional urban/indus trial pollution. Yoon et al. (2011b), using SeaWiFS 
data (with a different AOD retrieval algorithm) found positive trends in AOD during 
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springtime and fall over the Pearl River delta in China; positive trends in this region are also 
seen in this study, although only statistically-significant on an annual, rather than seasonal 
scale due to frequent cloud cover. In terms of overall trends, it is interesting to note that both 
East and South Asia exhibit similar patterns during the high aerosol loading seasons, i.e. a 
systematic increasing trend during the winter haze period and a rather variable AOD 
distribution during the dust-dominant spring/pre-monsoon season (Fig. 10). A recent study 
also indicates upward trends in fine-mode anthropogenic aerosols over China and India using 
inventory emissions data for sulfate and carbonaceous aerosols (Lu et al., 201 1). 

4.4 Zonal average trend 

In order to help study the radiative effects of aerosols on Earth’s climate, we also compute 
zonal-averaged trends (averaged from 180°W to 180°E over both land and ocean) based upon 
this 13 year SeaWiFS data record; these trends are displayed in Figure 12 (right panel). 
Trends from zonal mean of AOD anomalies for both annual and seasonal averages are 
calculated and shown from 55°S to 55°N at 5° intervals. The corresponding zonal averaged 
absolute AOD values are also depicted in Figure 12 (left panel) for comparison. As expected, 
the absolute AOD is much higher in the northern hemisphere compared to the southern 
hemisphere, which is mostly dominated by the oceans and is, in general, pristine. Most of the 
tropospheric aerosol burden resides between the equator and 40° N for all four seasons, 
except for the springtime Asian dust outbreak which peaks at mid-latitudes in the northern 
hemisphere, and the biomass burning activities over South America, Africa, and Southeast 
Asia in the southern hemisphere during the fall months (SON). 

As shown in Figure 12 (right panel), the corresponding overall annual trends are mostly 
positive over the latitude zone from equator to 30° N, where most of the aerosol sources 
reside, although there are significant self-compensating positive and negative trends within 
this zone as depicted in Figure 8. On a seasonal basis, the primary peak tendency in AOD 
occurs during the winter months around 25° - 40° N and is largely associated with changes in 
India and China as well as adjacent waters. However, the springtime zonal mean trends are 
slightly negative around 35° - 45° N due to weak decreases over the dust source regions in 
northwestern China, while remaining positive in the 10° - 30° N latitude zone as a result of 
larger increasing trends in the Arabian Peninsula dominating over the decreasing trend in 
Central America. The net positive summer trends over the tropics and subtropics in the 
northern hemisphere are also a product of the strong increasing Arabian trends and the weaker 
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decreasing trends in the Saharan dust dominant Atlantic Ocean. Compared to those in the 
northern hemisphere, the aerosol emissions in southern hemisphere, which are mostly related 
to biomass burning activities, peak in the fall (SON) months. Most of the aerosol tendency in 
this region is not statistically significant, especially near the smoke sources over land, due to 
large year-to-year fluctuation. However, a significant positive trend around 35° S is detected, 
which may be linked to the increasing amount of biomass burning smoke transported out of 
southern Australia and southeastern Africa. Overall, significant upward trends in annual mean 
AOD are found over the tropics and sub-tropics; they are largely amplified by the increasing 
dust emissions and outflow from the Arabian Peninsula during boreal spring and summer 
seasons. 

5 Conclusions 

Compared to other satellite sensors, SeaWiFS has unprecedentedly high accuracy and long- 
term stability (1-2% and 0.3%, respectively) in radiometric calibration as a result of rigorous 
lunar, solar and vicarious calibrations performed by the OBPG team. Therefore, the 
corresponding long-term drift in SeaWiFS AOD time series due to sensor calibration is 
negligible. However, as described in Section 1, several other factors in addition to calibration, 
such as retrieval algorithm deficiency, could also contribute to the temporal variability of 
satellite retrieved AOD data records. In order to examine the effect of retrieval algorithm on 
trend analysis, we compare the SeaWiFS v003 AOD products against ground-based 
AERONET measurements acquired over the 13 year mission lifetime. The results indicate 
that there is no apparent year-to-year shift over either land or ocean between these two 
measurements, suggesting that the SeaWiFS v003 data are suitable for estimating the 
interannual changes in AOD on both global and regional scales. 

Our trend analyses based upon the SeaWiFS data from 1998 to 2010 show that the global 
annual trend of AOD during this period, although weakly positive, is essentially negligible 
when compared to the magnitudes of contributions from other factors, including large-scale 
meteorological events such as ENSO and NAO. Our estimated SeaWiFS trend over global 
ocean is comparable to and has the same sign as the trend derived from the MODIS sensors 
aboard Terra and Aqua. However, it is in opposite sign to AVHRR data during the 
overlapping years. On regional scales, the estimated trends in this study suggest that the AOD 
tendency could be significantly large. We note that there are decreasing trends over the 
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eastern US and Europe, possibly due to a corresponding decrease manmade air pollution 
levels. Yet, over the rapidly developing countries such as China and India, significant 
increasing trends in AOD are seen in these source regions and their surrounding downwind 
oceans, particularly during the dry winter/post-monsoon months when the atmosphere is 
relatively stable, thus favoring accumulation of aerosols. 

For the mineral dust-dominated parts of the world, strong positive trends are detected over the 
Arabian Peninsula and the adjacent waters. In contrast, a negative tendency is observed in the 
emission and export of Saharan dust over the western North Africa and the North Atlantic. 
Overall, based on 13 years of data, a relatively contrasting pattern in trends seem to have 
emerged in the tropics/sub tropics largely modulated by dust emissions and transport processes 
encompassing the Saharan arid lands and the Arabian Peninsula and their downwind oceanic 
regions, with downward and upward tendencies, respectively. 

In order to reduce the uncertainty in estimating climate forcing caused by tropospheric 
aerosols, it is important to produce the best quality climatology of aerosol properties with 
sufficient data length by merging aerosol information obtained from the multiple satellite 
sensors that are currently available. SeaWiFS data extends from 1997 to 2010, which has 
many overlapping years with MODIS, MISR, AVF1RR, and TOMS, and thus provides a good 
benchmark to characterize instrument drifts due to calibration for other satellite 
measurements. However, 13 years of satellite measurements is not enough to establish long- 
term trends of aerosols, particularly in the tropics where the ENSO and other large-scale 
climatic factors could play a large role in modulating aerosol variability. Therefore, it is 
imperative to continue the high quality-EOS aerosol data records from SeaWiFS and MODIS 
with Visible and Infrared Imager Radiometer Suite (VIIRS) measurements onboard the 
currently operational Suomi NPP satellite and to be flown on future JPSS satellites to achieve 
the statistically adequate length of data records necessary for climate studies. 
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1 Table 1 . Comparisons of SeaWiFS-derived trends with those from AERONET measurements. 

2 ** indicates trends significant at 95%, * indicates trends significant at 90% 



Location 

Time 

Period 

AERONET 

SeaWiFS 

Trend 

(AOD/year) 

Std. 

Error 

Trend 

(AOD/year) 

Std. 

Error 

Ocean 

Capo Verde 

16N,22W 


-0.0040* 

0.0022 

-0.0040 

0.0031 

COVE 

36N,75W 


-0.0010 

0.0019 

-0.0026** 

0.0010 

Ascension 

Island 

7S,14W 


0.0012 

0.0013 

0.00083 

0.00073 

Naum 

0S,166E 


0.0013** 

0.00058 

0.00086 

0.00058 

Tahiti 

17S,149W 


0.00022 

0.00053 

0.00090** 

0.00019 

Shirahama 

33N,135E 


0.000063 

0.0023 

0.00056* 

0.00032 

La Parguera 

17N,67W 


0.0014 

0.00099 

0.00091 

0.00077 

Land 
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2 Figure 1. SeaWiFS average AOD at 550 nm from 1997-2010. 
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Figure 2. Yearly variations in the differences between SeaWiFS and AERONET AOD over 
(a) ocean and (b) land. Red is the mean difference in each year, while solid black is the 
median for each year. Dotted black are the 68% confidence interval. 
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2 

3 Figure 3. Correlation coefficients of SeaWiFS deseasonalized monthly AOD anomaly with 

4 the multivariate ENSO index (ME I) (a) from the same month, (b) leading by 6 months, and 

5 (c) leading by 1 year. Dots indicate significance at 95% confidence level. 
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2 Figure 4. Time series of AOD anomaly and ENSO Index over Indonesia region (80 S-80 N, 

3 900 -1300 E). 
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2 Figure 5. Lagged correlation of summertime average (JJA) mean AOD anomaly with winter 

3 (DJ) mean ENSO Index leading by 6 months. Dots indicate significance at 95% confidence 

4 level. 
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2 Figure 6. Correlation of AOD anomaly and NAO Index during winter (DJF) months. Dots 

3 indicate significance at 95% confidence level. 
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3 Figure 7. Comparisons of the absolute AOD (top) and AOD anomaly (bottom) averaged over 

4 global ocean as a function of year between SeaWiFS and various MODIS products. The 

5 estimated linear trend is 0.00098/year for SeaWiFS and 0.00061/year for MODIS Terra. 
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2 

3 Figure 8. Linear trend based upon deseasonalised monthly anomaly of AOD at 550 nm. Units 

4 are AOD/year. Dots indicate significance at 95% confidence level. The black and pink boxes 

5 and symbol (x) represent the ground based AERONET locations with AOD tendency of 

6 positive, negative, and negligible values, respectively. 


30 



DJF 


MAM 



3 Figure 9. Seasonal trends of SeaWiFS AOD anomaly from January 1998 to December 2010. 

4 Dots indicate significance at 95% confidence level. Units are AOD/year. The gray color 

5 represents the regions with insufficient sampling sizes for trend analysis due to the frequent 

6 coverage of clouds and snow/ice. 
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2 Figure 10. SeaWiFS time series of seasonal averaged AOD anomaly over northern India (20° 

3 - 30° N, 75° - 85° E), eastern China (30° - 40° N, 110° - 120° E), Arabian Peninsula (10° - 

4 35° N, 35° - 60° E), Europe (43° - 55° N, 0° - 30° E), and eastern United States (30° - 45° N, 

5 70° - 90° W). 
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2 Figure 11. Interannual variation of AOD (top panel) and AOD anomaly (bottom panel) from 

3 SeaWiFS and AERONET measurements co-located over the Solar Village site, for the period 

4 from February 1999 to April 2010. A systematic increasing aerosol loading signal is 

5 associated with the strengthening of the seasonal cycle, particularly during boreal spring and 

6 summer seasons, i.e. peak dust loading period. Interannual variations of Angstrom Exponent 

7 anomaly from AERONET data (bottom panel) also suggest a characteristic increase in coarse 

8 aerosol fraction. 
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2 Figure 12. Zonal average AOD trends from 55° S to 55° N (land and ocean). 


